Critical Coarsening without Surface Tension: the Voter Universality Class 
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We show that the two-dimensional voter model, usually considered to only be a marginal coarsen- 
ing system, represents a broad class of models for which phase-ordering takes place without surface 
tension. We argue that voter- like growth is generically observed at order- disorder nonequilibrium 
transitions solely driven by interfacial noise between dynamically symmetric absorbing states. 
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Coarsening phenomena occur in a large variety of situ- 
ations in and out of physics, ranging from the demixion of 
alloys [0 to population dynamics At a fundamental 
level, phase-ordering challenges our capacity to deal with 
nonequilibrium systems and our understanding of the 
mechanisms determining different universality classes. In 
many cases, phase competition is driven by surface ten- 
sion, leading to 'curvature-driven' growth. Coarsening 
patterns are then characterized by a single length scale 
L(t) ~ t x l z , where the exponent z only depends on gen- 
eral symmetry and conservation properties of the sys- 
tem Q. For instance, z = 2 for the common case of a 
non-conserved scalar order parameter (NCOP), a large 
class including the Ising model. In this context, the two- 
dimensional voter model (VM) |J, a caricatural process 
in which sites on a square lattice adopt the opinion of a 
randomly-chosen neighbor, stands out as an exception. 
Its coarsening process, which gives rise to patterns with 
clusters of all sizes between 1 and \ft (Fig. [l]) @||, is 
characterized by a slow, logarithmic decay of the den- 
sity of interfaces p ~ 1/lnt (as opposed to the algebraic 
decay p ~ 1/L(t) ~ t~ x l z of curvature-driven growth). 
The marginality of the VM is usually attributed to the 
exceptional character of its analytic properties [||-|| . 

In this Letter, we show that, in fact, large classes 
of models exhibit the same type of domain growth as 
the simple VM, without being endowed with any of its 
peculiar symmetry and integr ability properties. We ar- 
gue that voter-like coarsening is best defined by the ab- 
sence of surface tension and that it is generically observed 
at the transitions between disordered and fully-ordered 
phases in the absence of bulk fluctuations, when these 
nonequilibrium transitions are driven by interfacial noise 
only. Finally, we discuss the universality of the scaling 
properties associated with voter-like critical points. 

We first review the properties of the usual two-state 
VM, emphasizing those of importance for defining gen- 
eralized models. A 'voter' (or spin) residing on site x 
of a hypercubic lattice can have two different opinions 
s x = ±1. In any space dimension d, an elementary move 



consists in randomly choosing one site and assigning to it 
the opinion of one of its randomly chosen nearest neigh- 
bors (n.n.). This ensures that the two homogeneous con- 
figurations (where all spins are either + or — ) are ab- 
sorbing states, and that the model is ^-symmetric un- 
der global inversion (s x — > — s x ). Recasting the dynamic 
rule in terms of n.n. pair updating, a pair of opposite 
spins H — is randomly selected, and evolves to a ++ or 
pair with equal (1/2) probabilities. The rates of cre- 
ation of + and — spins being equal, any initial value of 
the global magnetization m is conserved in the limit of 
large system sizes. 

Another prominent feature of the VM is a 'duality' || 
with a system of coalescing random walks: going back- 
ward in time, the successive ancestors of a given spin 
follow the trail of a simple random walk (RW); compar- 
ing the values of several spins shows that their associated 
RWs necessarily merge upon encounter. This correspon- 
dence allows to solve many aspects of the kinetics of the 
VM, because it implies that the correlation functions be- 
tween an arbitrary number of spins form a closed hierar- 
chy of diffusion equations ||J^]. In particular, the calcula- 
tion of the density of interfaces p m (t) (i.e. the fraction of 
H — n.n. pairs) starting from random initial conditions 
(r.i.c.) of magnetization m, is ultimately given by the 
probability that a RW initially at unit distance from the 
origin, has not yet reached it at time t. Therefore, ow- 
ing to the recurrence properties of RWs, the VM shows 
coarsening for d < 2 (i.e. p m (t) — > when t — > oo). For 
the 'marginal' case d = 2 — on which we mainly focus 
henceforth — , one finds the slow logarithmic decay |5]||: 



Pm(t) = (1-m 2 ) 
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D being the diffusion constant of the underlying RW 
{D = 1/4 for the standard case of n.n., square lattice 
walks, when each spin is updated on average once per 
unit of time). This peculiar behavior, which contrasts 
with the algebraic decay encountered in NCOP growth, 
thus seems tantamount to the fact that d = 2 is the crit- 
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ical dimension for the recurrence of RWs. 

Beyond these analytic properties, the nature of the 
coarsening displayed by the VM in two dimensions is 
highlighted by studying the evolution of a 'bubble' of 
one phase embedded into another. For curvature-driven 
NCOP growth, the volume of a bubble decreases linearly 
with time under the effect of surface tension. In the VM, 
however, large bubbles do not shrink but slowly disinte- 
grate as their boundary roughens diffusively to reach a 
typical width of the order of their initial radius ro (Fig. [l|, 
top). As long as ^/t <C ro, conservation of the magnetiza- 
tion is still effective, and implies that radially- averaged 
magnetization profiles m(r, t) have a stationary middle 
point. This indicates that curvature plays no role, a 
fact further confirmed by the observation that the m(r, t) 
curves are the same as the transverse profiles of an infi- 
nite, straight interface, whose derivatives can be shown 
to be given by a simple Gaussian of variance 2Dt. 




FIG. 1. Domain growth in the usual VM (system size 
256 2 ). Top: snapshots at t — 4,16,64,256 during the evo- 
lution of a bubble of initial radius ro = 180 (thin circle). 
Bottom: same from symmetric r.i.c. 

We interpret the above behavior as characteristic of the 
absence of surface tension, a 'physical' property which, 
together with the fact that coarsening does occur (i.e., 
lim^oo p — 0), we conjecture to be constitutive of d = 2 
voter-like domain growth (as defined by the logarithmic 
decay of p m (t)). To investigate which of the properties 
outlined above (existence of two absorbing states, Z 2 - 
symmetry, duality with coalescing RWs, m-conservation) 
is necessary to produce voter-like growth, we now embed 
the simple VM studied above into large families of rules. 

Consider first that a H — pair is now only updated 

with probability p (Pr(H ► H — ) = 1 — p), and that 

Pr(± > ++) = pq and Pr(± > ) = p(l - q), 

so that for q ^ - the m-conservation is broken. Sup- 
pose next that p and q depend on the local configuration 
around the H — pair, for instance via n + and n~, the 
numbers of + (resp. — ) n.n. of the — (resp. +) site 
in the central pair. (Note that as soon as the transi- 
tion probabilities do vary with the local configuration, 
the duality with coalescing RWs breaks down.) Conser- 
vation of m (Pr(H ► ++) = Pr(H > )) is then 

simply expressed by the conditions q n +, n - = \, to be 



obeyed Vn + , n G {1,2,3,4}, whereas p n - }7l + = p n +, n - 
and q n -, n + = 1 — q n +,n- stand for ^-symmetry. 

With the exception of the usual VM (p = 1, q — 1/2), 
the m-conserving and ^-symmetric rules, defined by 
Qn+,n- = 1/2 andp n - n + =_p n +,n-? do not benefit from 
the duality property used above to derive, e.g., Eq. ([!]). 
Nevertheless, we have found (by coarsening, bubble and 
line experiments ||D that they behave like the simple 
VM, but with a diffusion constant D ^ |. Thus, integra- 
bility is not a necessary condition for voter-like growth. 

The subset of m-conserving rules without Z 2 - 
symmetry is of particular interest: no surface tension 
is expected, but one may wonder about the effects of 
the asymmetry on growth properties. Below, we use 
the case p n +, n - = l/n + , which is nothing but the in- 
terface growth model introduced by Kaya, Kabakcioglu, 
and Erzan (KKE) at its 'derealization' transition [fzj. 
The asymmetry of the rule manifests itself at the mi- 
croscopic level (by looking at local 3-spins configurations 
||), but also at the macroscopic level, by the evolution 
of the magnetization profiles in 'line' or 'bubble' experi- 
ments. As for the usual VM, the profiles obtained in both 
cases are identical, their width scale as 2DkkeV^ with 
Dkke = 0.36(2) ||, but they are asymmetric (Fig. ^|a). 
Voter-like behavior is further confirmed by the study of 
phase ordering following r.i.c. of magnetization m. The 
density p m (t) displays the signature scaling spelled out 
by Eq. (g) (Fig. |b), and the ensuing value D = 0.34(1) is 
in fair agreement with the value Dkke determined above, 
which suggests that the KKE rule asymptotically behaves 
like the usual VM with an effective diffusion constant. 
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2. Domain growth in the KKE model and the n.n. 
(D = |) VM. (a): derivative of rescaled magnetization pro- 
files at times t = 16, 64, 256 (data from a bubble experiment; 
initial conditions: m(r, 0) = sign(r — ro) with ro = 4096). 
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Dashed line: exact result for the VM dm/dy = (2/y / 7r)e _y 
with v—{r — ro)/y/i; (b): scaling of p m {t) (as suggested by 
Eq. (pi)) in coarsening experiments for m — 0, ±0.2, ±0.4 in 
systems of size 16384 2 ; dashed line: exact result for the VM. 



All the other m-conserving, asymmetric rules that we 
have studied reveal a similar voter-like behavior. We con- 
jecture that they almost all do so, as long as the p n +, n - 
are all strictly positive (to avoid possible 'blocking' in 
certain configurations). In other words, m-conservation 
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appears as a sufficient condition for voter-type growth, 
even in the absence of ^-symmetry or integrability. We 
now show that m-conservation is not a necessary condi- 
tion since some non-conserving ^-symmetric rules with 
two absorbing states exhibit voter-like features. 

If pair-update rules are convenient to control m- 
conservation, single-site update models suffice to study 
rules which need not possess this property. We thus con- 
sider a family of 'kinetic Ising models', in which a spin 
is flipped with a probability r 3i h depending on its value 
s = ±1 and on its local field h G {-4,-2,0,2,4}, i.e. 
the sum of the values of its four n.n.. The ^-symmetric 
rules of this family, obeying r- s -h = r 3j h, can be param- 
eterized by five transition probabilities ru = ri 5 ^. Such 
rules are usually not integrable. Again, the VM, given by 
rh = \ — f, is an exception. Its marginal situation in this 
context has already been noted [jTo| ? |lTj| in the particular 
case where the transition probabilities only depend on h, 
a condition expressed here by r- S: h = 1 — r s ,h- One is 
then left with only two free parameters, p\> = pi = r<i, 
which can be respectively interpreted Q as a measure 
of the strength of 'bulk' and 'interfacial' noise. In the 
(PbPb) plane, the VM {p\ = \, Ph = 0) turns out to be 
the endpoint of a line of Ising-like (second-order) phase 
transitions |n],|ll| , at the border between the ordered and 
disordered phases in the absence of bulk noise (pb = 0). 

We claim that this situation is generic: the 4- 
parameter space of the general kinetic Ising model above 
without bulk noise (space defined by = 0) is divided 
into a disordered and an ordered region separated by a 
codimension-1 critical manifold of voter-like rules |T^ ]. In 
Fig. ||, we show a section of this manifold in a plane which 
contains the VM. Its location was determined by study- 
ing the evolution of p following r.i.c. with zero magneti- 
zation at various parameter values. Crossing the critical 
manifold, the characteristic logarithmic decay p oc 1/lnt 
separates the region where p oc 1/y/t (curvature-driven 
growth typical of the ordered phase) from the disordered 
region where p saturates to finite values (Fig. ||b). A pe- 
culiar feature of this voter-like critical manifold is that, 
using the language of the renormalization group (RG), 
m = is a 'weakly attracting' fixed point of the dynam- 
ics: starting from any mo ^ 0, ±1 (mo = m(0)), one ob- 
serves that m(t) cx 1/lnt (Fig. ||c). Yet, the decay of the 
interface density behaves like in the usual VM (Fig. ||d), 
even for mo ^ 0, because the slow evolution of m(t) intro- 
duces only 0(1/ In 2 t) corrections to Eq. (|l|) and therefore 
does not alter the value of the effective diffusion constant. 
This also explains why our bubble and line experiments 
— which strictly speaking probe only a zero-measure set 
of random initial conditions — performed for rules on the 
voter-like manifold, do give the same results as for the 
usual VM ||, in spite of the absence of m-conservation. 

Similar investigations have been undertaken in various 
planes of the 4-parameter space, always yielding the same 
results. This suggests to conjecture that all critical Z2- 



symmetric rules without bulk noise form a codimension-1 
'voter-like' manifold separating order from disorder, char- 
acterized by the logarithmic decay of both p and m. 

Does the above picture extend to non-conserving, 
asymmetric rules without bulk noise? In this case, even 
though two absorbing states exist, their dynamical roles 
are not symmetric, and such order-disorder transitions 
are known Jl3|,|l4| to belong generically to the univer- 
sality class of directed percolation (DP). We have indeed 
observed |9| that for such rules both m and p scale at crit- 
ically with the same exponent (3 BP /VP P = 0.450(1) [El. 



(a) 



0.5 



1 






A - 




z 2 VM 






1 


1 1 





2 4 6 2 4 

FIG. 3. Voter-like rules not conserving m: (a): phase di- 
agram in the plane (r_2,r_4) for ro = \,r<i — r-4/4. (b): 
pit) from mo — r.i.c. increasing r_4 (from top to bottom) 
with r-2 — 0.275 (dashed line in (a)). Inset: interface den- 
sity voter- like decay at the critical point (system size 2048 2 ). 
(c,d): m(t) (c) and p(t) (d) from r.i.c. with various mo at the 
critical point described in (b) (system size 16384 2 ). 



Having specified the required conditions to observe 
voter-like coarsening, we now turn to the investigation of 
the scaling laws associated to voter-like points and to an 
assessment of their universality. To our knowledge, there 
has been only one attempt to determining the scaling ex- 
ponents governing the approach to the usual VM. Within 
the reduced (pi,pt>) parameter plane defined above, the 
authors of [[10| noted first that increasing p[ along the 
zero-bulk- noise line p\> = 0, up to the VM point = |), 
the (static) magnetization m s jumps from to ±1. Thus, 
the order parameter exponent f3 = 0, and, for such a 
first-order d = 2 nonequilibrium critical point, the sus- 
ceptibility and the correlation length exponents should 
satisfy 7 = 2v. A finite-size scaling (FSS) analysis of the 
fluctuations of m s as £ = £>i — \ ^0~ confirmed the lat- 
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ter relation, and gave 7 ~ 1.25. However, repeating these 
simulations with much better statistics, we evidenced ||] 
a systematic decrease of the local exponent 7 as £ — > 0~. 
On general grounds, corrections to scaling in this problem 
are expected to be logarithmic [|l5| and, thus, the above 
7 value must be taken cautiously. In fact, a more reliable 
approach of the VM is from the ordered side (e — > + ), 
along the zero-bulk noise line. The curvature-driven- 
growth regimes p ~ £(e)/L(t) ~ e~ v /t 1 / 2 , which eventu- 
ally settle then in large enough systems, are less prone to 
the logarithmic corrections mentioned above. The scaling 
of the correlation length £(e) yields v = 0.45(7), compat- 
ible with the simple 'diffusive' value v = \ (Fig. [|). 
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FIG. 4. Phase-ordering from m — r.i.c. approaching the 

VM from the ordered side (system size 8192 2 ). (a) pvsVt 

for various e between 10 -3 and 5.10 -2 (from bottom to top); 

(b) extracted correlation length £ vs y/e. 

Compounding the above results, a conservative in- 
terpretation of the data is that (3 = 0, 7 = 1, and 
v = ^, with appropriate logarithmic corrections to scal- 
ing. This is further confirmed by similar, preliminary, 
results obtained when approaching any of the voter-like 
models studied above ||. We finally note that these 
exponents would obey all the standard scaling relations 
valid at a critical point, such as (3 = [d — X(H)]/X(Ti), 
7 = [2X(H) - d]/X(Ti), or v = l/A^), with (in d = 2) 
X(Ti) = X(H) = 2. We tentatively associate the for- 
mer RG-eigenvalue with the relevance of inter facial noise, 
and the latter with the presence of a 'dynamically self- 
induced' magnetic field breaking the ^-symmetry (as 
happens in KKE-like models). The behavior of the 
classes of voter-like models we have defined is consistent 
with such a set of exponents, which would thus charac- 
terize (in d = 2) the genuine voter critical point. For 
instance, converting the static critical behavior m ~ 
to a dynamical one via £ ~ e~ v ~ L(t) ~ t x l z gives m ~ 
t -P/vz m Now? if p = strictly, this implies m(t) = Cst, 
while if (3 = + is interpreted, as is customary, including 
logarithmic corrections, then m(t) ~ 1/lnt. It is also 
nice to note that, despite the presence of the a priori rel- 
evant eigenvalue A (if) = 2, m-conserving rules without 
^-symmetry may still display voter-like growth, because 
the critical exponent 5 (which describes how the order 
parameter behaves under a magnetic field: m ~ H 1 ' 5 ), 



would be formally given by 1/5 = —1 + d/X(H) = 0. 

Naturally, the RG picture sketched above needs to be 
substantiated by the study of an appropriate field the- 
ory, an endeavor left for future work. In particular, one 
would like to have a better understanding of the RG flow 
around the voter critical point as the space dimension 
varies, and of the special role played by d — 2, which 
appears as the upper critical dimension where voter-like 
phase-ordering can be induced by the sole presence of in- 
ter facial noise [[l6|]9|]. At any rate, we believe the results 
presented here reveal that the voter model, which was, up 
to now, perceived as a marginal system, embodies in fact 
a broad class of models for which coarsening occurs in the 
absence of surface tension. Our findings also suggest that 
the critical behavior associated to this class of systems 
is perhaps best characterized as an order-disorder transi- 
tion driven by the interfacial noise between two absorbing 
states possessing equivalent 'dynamical roles', this sym- 
metry being enforced either by ^-symmetry of the local 
rules, or by the global conservation of magnetization. 
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